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I.  INTRODUCTION 


This  final  report  contains  the  results  of  the  project  SPC-99-4060  and  concentrates  mainly 
on  a  summary  of  high  temperature  defect  structure  of  CdTe  and  defect  modeling.  Some  re¬ 
sults  already  presented  in  the  Progress  Report  submitted  in  January  2000  were  extended 
and  an  improved  evaluation  and  discussion  of  experimental  data  was  performed.  High  tem¬ 
perature  defect  equilibrium  was  studied  on  both  the  Cd  and  Te  sides  of  the  P-T  diagram 
and  a  complex  theoretical  approach  is  presented.  The  main  effort  concentrated  on  iu-situ 
galvanomagnetic  measurements  in  the  liquid  CdTe,  where  in  the  case  of  measurements  of 
Hall  effect  no  experimental  data  were  published  so  far.  The  results  of  crystal  growth,  char¬ 
acterization  of  crystals  grown  under  different  conditions  by  optical  methods  and  modeling 
of  crystal  growth  were  given  in  the  Progress  report. 

The  wide-bandgap  binary  semiconductor  CdTe  (CT)  with  the  bandgap  1.5  eV  at  room 
temperature  is  used  as  a  source  material  for  preparation  of  detectors  of  gamma  radiation, 
as  a  substrate  for  fabrication  of  detectors  of  infrared  radiation  on  the  basis  of  (CdHg)Te 
and  also  as  a  material  for  electrooptical  modulation  of  signals.  The  recent  progress  in  high 
technology  brings  a  renewed  interest  to  the  research  of  basic  properties  of  CT  [1,2]. 

In  order  to  optimize  the  function  of  the  devices  it  is  necessary  to  improve  both  structural 
and  material  parameters  of  the  single  crystalline  CT.  According  to  the  state-of-the-art  of 
technology  of  growth  of  CT  single  crystals,  both  demands  lead  to  minimization  of  influence 
of  native  point  defects,  which  can  affect  e.g.  the  room  temperature  mobility  and  lifetime 
of  carriers  or  the  concentration  of  precipitates  of  one  of  the  components  (Cd,Te)  at  higher 
defect  concentrations.  To  prepare  the  high  quality  single  crystalline  CT  it  is  necessary  to 
know  the  equilibrium  concentrations  of  point  defects  and  diffusion  processes,  which  establish 
the  equilibrium  concentrations  during  cooling/heating  or  when  the  pressure  of  one  of  the 
components  above  the  crystal  is  changed.  The  real-time-diffusion  occurs  at  temperatures 
above  500°C. 

The  equilibrium  concentration  of  electrically  active  point  defects  in  the  whole  region  of 


1 


stability  of  the  solid  phase  depends  on  the  temperature  and  pressure  of  one  of  the  com¬ 
ponents.  It  implies,  that  the  concentration  of  free  carriers  in  the  sample  changes  in  the 
dependence  on  both  the  temperature  and  the  pressure  of  one  of  the  components,  e.g.  Cd, 
as  well.  The  equilibrium  concentrations  of  defects  can  be  experimentally  studied  in-situ  at 
high  temperature  by  the  measurement  of  the  Hall  coefficient  7?//  and  conductivity  o ,  from 
which  the  Hall  mobility  ////  =  Rhc t  is  calculated. 

The  measurement  of  Ru  and  a  can  be  performed  by  two  ways.  In  the  first  way  the 
sublimation  of  the  more  volatile  component  (Cd)  is  prevented  by  the  sample  encapsulation, 
e.g.  by  SHN4  deposited  in  a  pyrolytic  reactor.  The  content  of  Cd  and  Te  does  not  change 
and  only  the  temperature  dependence  is  studied  in  this  case.  In  the  second  method,  which 
is  used  in  our  case,  the  content  of  Cd  in  the  sample  changes  due  to  the  sublimation  from 
the  free  surface.  The  equilibrium  is  controlled  here  by  a  proper  Cd  overpressure  which  is 
established  in  the  two-zone  furnace  by  a  temperature  of  the  Cd  source.  The  advantage  of 
this  arrangement  is  the  possibility  to  measure  properties  of  CT  in  the  phase  P-T  diagram 
both  in  intrinsic  conditions  and  under  Cd  or  Te  overpressure.  The  measurement  of  II//  and 
a  is  performed  in  dependence  on  the  pressure  and  temperature  of  Cd  when  the  equilibrium 
concentration  of  electrically  active  defects  is  established  by  diffusion  processes.  Experimental 
data  show  as  expected,  that  R//  and  a  depend  both  on  pressure  of  Cd  and  temperature,  but 
the  [in  is  near  intrinsic  conditions  practically  independent  of  Cd  pressure  P It  implies 
that  there  is  not  important  impurity  scattering  in  this  case.  Figure  1  shows  the  phase 
P-T  diagram  for  CT  with  a  marked  region  of  stability  of  CT,  where  our  measurements  were 
performed  (Cd-rich  area  of  the  phase  diagram). 

II.  EXPERIMENT 

The  measurement  was  performed  on  samples  from  single  crystals  prepared  by  vertical 
cooling  in  the  temperature  gradient.  The  high  temperature  measurements  were  done  in  a 
quartz  ampoule  with  a  diameter  17  mm  located  in  a  furnace 
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Pig.  1 .  The  P-T  phase  diagram  of  CdTe  with  the  marked  area  (dotted)  at  which  the 
measurements  were  performed.  Ps,  P'7,  and  PL  plot  the  solidus,  gasus,  and  liquidus  stoi¬ 
chiometry  condition,  respectively.  P°  is  the  Cd  partial  pressure  over  pure  Cd. 
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Fig. 2.  Experimental  and  theoretical  data  of  Hall  mobilities  in  CdTe.  The  low  temper¬ 
ature  points  (o)  are  according  to  Ref.  [3].  Our  high  temperature  experiment  (•)  is  shown 
together  with  previous  measurements  (A)  [6]  and  (+)  [7].  Our  data  confirm  the  measure¬ 
ments  of  [7].  The  full  line  draws  the  best  theoretical  fit.  The  dash  line  is  without  effect  of 
the  L- valley.  The  dotted  line  is  for  temperature  independent  A E  =  0.19  eV.  The  inset  plots 
the  scheme  of  the  model  band  structure. 
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Fig. 3.  The  Hall  mobilities  within  different  types  of  scattering  included  (left  axis,  full) 
and  the  conduction  band  Hall  factors  (right  axis,  dash).  The  mobility  pr  without  effect  of 
the  L-valley  -  (line  a).  The  mobility  /ip  including  the  intervalley  scattering  -  (line  b).  The 
total  Hall  electron  mobility  -  (line  c).  The  mobility  -  (line  d).  rr  -  (line  r),  tl  -  (line  L). 
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It  was  possible  to  keep  the  sample  temperature  with  the  precision  of  about  1°C  and  t  o 
change  the  gas  pressure  of  one  of  the  components  (in  our  case  Cd)  in  the  interval  2  x  10  4 
2  atm.  The  sample  with  tungsten  or  molybdenum  spring  contacts  in  the  van  der  Pauw 
configuration  was  placed  in  a  quartz  holder  together  with  the  Pt-PtRh  thermocouple. 

III.  HIGH  TEMPERATURE  MOBILITY 

A.  Model  and  Theory 

We  shall  concentrate  on  electrons  in  this  report,  because  the  influence  of  holes  is  weak 
and  can  be  included  in  a  simple,  standard  way.  The  basic  problem  at  the  evaluation  of  the 
drift  electron  mobility  /i  at  temperatures  (500-1000°C)  is  the  determination  of  scattering 
mechanisms,  which  participate  in  scattering  of  current  carriers.  A  number  of  detailed  exper¬ 
imental  measurements  of  mobility  was  performed  at  low  temperatures  (50-30QK),  which  can 
be  explained  well  by  scattering  of  electrons  on  ionized  impurities  and  polar  optical  phonons 
[3].  The  influence  of  ionized  impurities  on  scattering  above  400  K  is  negligible  in  pure  sam¬ 
ples  and  the  dominating  scattering  on  longitudinal  optical  phonons  is  expected  describing 
/x  in  the  form 


where  the  effective  mass  rn  =  0.096m0,  the  longitudinal  optical  phonon  energy  fiui0  — 
21.0  meV,  and  the  optical  and  static  dielectric  constants  eopt  —  7.1  and  es  —  10.3  are 
material  parameters.  The  dimensionless  function  G  is  calculated  according  to  Ref.  [4]  and 
K  describes  the  correction  to  nonparabolic  conduction  band  within  the  Kane  model  [5].  The 
mobility  is  only  slightly  dependent  on  uq  at  high  T. 

The  data  of  measurements  at  high  temperatures  (T  >  500°C)  published  so  far  [0,7] 
exhibit,  however,  substantially  stronger  decrease  of  /i//  which  reaches  above  900°C  less  than 
one  half  of  its  theoretical  prediction  [8].  Furthermore  values  from  Ref.  [7]  are  «  20%  lower 
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than  those  of  Ref.  [6].  The  verification  of  the  experimental  data  and  the  explanations  of  t  he 
effect  is  the  main  topic  of  this  report. 

Onr  measurements  of  Hall  coefficient  R//  and  electric  conductivity  a  confirm  the  relia¬ 
bility  of  values  published  in  Ref.  [7]  (see  Fig.  2). 

There  are  two  recommendations  found  in  the  literature  to  explain  the  fast  decrease  of 
fin.  The  first  one  [6]  expects  a  strong  temperature  dependence  of  material  parameters, 
particularly  the  static  dielectric  constant.  In  agreement  with  Ref.  [8]  we  do  not  approve 
such  explanation.  In  view  of  (1)  and  the  values  of  the  material  parameters,  es  would  have 
to  increase  twice  within  the  temperature  interval  300°C  -  1000°C.  Such  behaviour  of  es  is 
not  observed  in  semiconductors.  Also  in  case  of  Copt  and  in  their  temperature  evolution  act 
probably  more  to  the  higher  than  to  the  lower  /x//. 

The  second  model  explanation  is  based  on  an  assumption  of  multivalley  conduction 
which  comprises  both  rlc  minimum  and  four  Llc  minima  of  the  Brillouine  zone  [9].  The 
heavier  electrons  in  the  L-minima  reduce  /x//  comparing  to  /xp  of  pure  T-valley  electrons.  The 
similarity  to  the  same  effect  in  GaAs  suggests  that  the  T \c  to  L\c  separation  AE  may  be  much 
less  than  ~1  eV  which  is  obtained  by  band  structure  calculations  [10,11].  Though  the  idea  is 
known  for  a  long  time,  there  is  no  concrete  calculation  of  the  high-teinperature  transport  in 
CT  except  the  numerical  modeling  of  the  Gunn  effect  and  hot-electron  phenomena  [12,13]. 
We  mean  in  this  report  to  complete  the  knowledge  about  CT  in  this  respect  calculating  the 
electron  mobility  in  the  rlc-valley  /xp  and  Llc-valleys  /xL,  as  well  as  the  total  Hall  mobility  /x// 
at  high  temperature  including  the  Fic+Llc  intervalley  scattering  and  multivalley  conduction. 

Our  treatment  is  based  on  the  solution  of  the  Boltzmann  transport  equation  within 
the  relaxation  time  approximation.  The  polar  optical  phonon  (PO)  and  acoustic  phonon 
(AC)  as  dominant  intravalley  and  intervalley  scatterings  [12]  are  used.  Due  to  the  intention 
to  study  transport  at  high  T,  this  approach  can  be  applied  also  for  the  optical  phonon 
scattering.  To  reduce  the  amount  of  unknown  parameters  all  valleys  are  assumed  spherical. 
The  corresponding  relaxation  times  Tpo  for  polar  scattering  and  tac,  for  accoustic  scattering 
read 
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(2) 


,  _  27re0  £s£opt  exP  ftr  “  1  ^  (  huo  \ 

rpo(E)  -  #*£*.-*«  ^  ' 


tac(E)  = 


hc-i 


(3) 


7T  E%kBTD{E)' 

whore  Ak  expresses  the  single  valley  dispersion,  q  =  7  x  1010  Nm-2  is  the  longitudinal  elastic 
constant  [14],  and  D  =  Dp  +  4 DL  is  the  total  conduction  band  density  of  states  formed  by 
the  single  T  and  four  L  valleys,  as  well.  The  intervalley  deformation  potential  is  not  known 
in  CT,  thus  we  use  the  intravalley  value  E\  =  4eV  [14].  This  is  also  the  reason  why  we  do 
not  involve  other  types  of  scattering  into  the  calculations.  Comparing  to  the  uncertainty  of 
the  used  material  parameters  such  an  effort  would  not  provide  credible  improvement  of  t  he 
results.  The  function  G  in  (2)  expresses  the  correction  to  inelastic  PO  scattering.  It  allows 
to  use  Tpo  also  at  lower  temperature  where  kBT  ~  Fiui0.  However,  such  Tpo  produces  at  low 
T  only  approximate  Hall  factor  rB.  The  intravalley  Tpo  can  be  completed  by  the  intervalley 
PO  scattering  terms 


W* )  =  4 5TP) 


tpo,l(E)  — 

td  = 


2  td 


2Dr(E)  +  3DL,(E) 

3TT£0h  ^s^opt 


(4) 

(3) 


4e2kBTa2 es  -  eopt' 

which  yield  the  relaxation  time  comparable  with  tac ■  The  k-space  intervalley  separation 
is  involved  through  lattice  constant  a  =  6.48A.  The  total  relaxation  time  for  all  valleys  is 
given  as 


--Ef1. 


(6) 


where  T{  label  all  relevant  scatterings  (2)-(5)  which  yields  the  T  and  L  mobilities.  The  Hall 
mobility  /i//  is  then  obtained  in  the  form 

rr/ifor  +  rL\i]nL  -  rpfi2pp 


Hh 


/ipnr  +  ^LnL  +  HPP 


(7) 
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whore  rp,  r^,  and  rp  are  the  Hall  factors  <  r2  >  /  <  r  >2  and  nr,  ni,  and  p  are  respective 
electron  and  hole  concentrations,  which  are  obtained  in  an  obvious  way.  We  found  that  the 
fit  quality  is  not  influenced  significantly  by  a  choice  of  vnL.  We  use  mL  =  0.35m0  in  this 
report.  This  value  was  successfully  used  in  numerical  simulations  of  the  Gunn  effect  [12], 

B.  Results  and  discussion 

The  valley  separation  A E  is  fit  to  the  experimental  data.  We  analyze  temperature 
dependent  A E  =  AE0  +  aT  which  results  to  the  best  fit  in  the  form 

AE  =  0.29  -  10-4T(A')  (eV),  (8) 

The  results  of  the  calculations  together  with  experimental  data  are  presented  in  Fig.  2.  We 
see  that  eq.  (8)  describes  experimental  data  quite  well.  For  an  estimation  of  the  influence 
of  AE  to  p//,  the  theoretical  data  for  a  constant  A£'=0.19  eV  are  shown  by  the  dotted  line 
in  Fig.  2.  If  we  consider  possible  deviations  of  material  parameters  and  experimental  errors, 
we  estimate  the  precision  of  our  fit  of  AEo  to  ±0.05  eV.  The  T- independent  AE  produces 
worse  fit  with  a  deviation  from  the  experiment  about  three  times  greater  than  the  fit  with 
T-dependent  AE.  Within  the  estimated  experimental  precision  we  cannot,  however,  confirm 
a  <  0  definitely.  The  a  =  0  produces  the  fit  which  expresses  the  basic  character  of  Hn(T)  as 
well.  On  the  other  way,  a:  >  0  results  in  the  fit  with  significant  deviations  from  experiment. 

For  a  comparison,  /i//  with  L- valley  excluded  is  plotted  in  Fig.  2,  too.  In  Ref.  [7],  the 
drift  mobility  from  the  Hall  mobility  using  acoustic-phonon  Hall  factor  3n/8  was  obtained. 
In  view  of  our  Hall  mobility  presentation  we  multiplied  the  data  of  [7]  by  that  value  to  return 
to  ph  here.  Numerical  simulation  of  the  Gunn  effect  [12]  were  performed  with  AJF=0.51 
eV  and  mi  =  0.35mo-  This  AE  is  however  too  high  to  explain  our  transport  data.  For  a 
comparison  A£'=0.3  eV  and  mi  =  0.22mo  are  used  in  GaAs.  The  value  of  a  in  eq.(8)  looks 
reasonable  in  respect  to  the  energy  gap  [15] 

Ea  =  1.G22  -  3.5  x  1(T4T  -  1.1  x  1(T7T2  {eV).  (9) 
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The  related  ellipsometry  measurements  and  theoretical  calculations  of  a  due  to  thermal 
expansion  and  life-time  broadening  were  reported  for  GaAs  [16].  The  results  show  that  a 
should  be  negative  in  GaAs  as  well,  especially  due  to  the  life-time  broadening  effect. 

The  effect  of  the  multivalley  conduction  and  the  intervalley  scattering  is  demonstrated 
in  Fig.  3.  We  can  distinguish  two  effects,  which  reduce  comparing  to  the  single  T-valley 
transport  represented  by  the  line  (a).  The  electrons  in  the  T-valley  at  the  energy  above  A E 
scatter  strongly  into  L-valleys  due  to  the  high  density  of  states  there.  If  these  high  energy 
levels  are  occupied,  the  mobility  /zp  is  depressed.  Such  a  phenomenon  is  represented  by  the 
line  (b).  The  heavy  electrons  in  the  L- valley,  shown  by  line  (d),  are  less  mobile  but  their 
concentration  quickly  increases  at  high  T.  Finally,  due  to  the  stronger  effect  of  transport  in 
the  L- valley  to  the  denominator  than  to  numerator  in  (7),  the  total  Hall  mobility  drops  to 
its  final  shape  as  the  line  (c). 

The  involvement  of  the  Hall  factors  in  (7)  influences  /z//  significantly.  As  it  is  shown  in 
Fig.  3,  due  to  the  intervalley  scattering  the  relaxation  time  is  reduced  at  energy  above  A E 
which  results  in  the  increased  rp. 


IV.  HIGH  TEMPERATURE  DEFECT  STRUCTURE  -  SOLID  STATE 

A  number  of  galvanomagnetic  (GM)  measurements  performed  on  samples  quenched  after 
annealing  under  defined  pressure  of  one  of  the  components  [15]  or  in-situ  at  high  tempera¬ 
tures  (HT)  [6,7,17  19,22,24]  exists,  which  were  used  to  establish  equilibrium  concentrations 
of  electrically  active  defects  in  the  region  of  stability  of  solid  phase.  The  former  method 
cannot  guarantee  ideal  quenching.  Due  to  a  low  emission  of  radiation  at  lower  tempera¬ 
tures  and  a  low  heat  conductivity  of  CdTe  only  a  surface  layer  of  samples  is  cooled  quickly 
and  a  concentration  gradient  from  surface  to  the  middle  of  the  sample  is  formed.  At  high 
temperatures  the  emission  of  radiation  and  therefore  also  the  cooling  is  quicker,  but  as  a 
result  of  the  high  diffusion  coefficient  of  Cd  the  quenched  room  temperature  concentration 
can  correspond  to  a  unknown  lower  than  expected  annealing  temperature.  In  the  case  of 
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all  in-situ  measurements  published  so  far  the  electroneutrality  condition  (ENC)  which  wras 
used  for  high  temperature  (HT)  equilibrium  of  defects  and  the  interpretation  of  obtained 
values  of  the  Hall  coefficient  R h  and  conductivity  a  were,  however  not  correct. 

We  calculate  the  HT  electron  transport  properties  for  CdTe  with  respect  to  the  combined 
transport  both  in  Tic  and  in  L\c  space.  The  nonparabolic  kp  band  approximation  within  the 
Kane  model  is  used  for  Tlc  minimum  with  the  effective  mass  mr  =  0.09Grao  [30].  The  band 
structure  near  the  Llc  point  is  not  known  precisely.  We  use  parabolic  band  with  effective 
mass  mi  —  0.35mo  here  and  the  T\c-L\c  valley  separation  A E  =  0.29  —  10~4T(R')(eV), 
which  was  obtained  from  the  fit  of  experimental  high  temperature  mobility  [57]. 

Calculation  of  hole  transport  properties  is  more  simple,  when  mostly  heavy  holes  with 
effective  mass  in/,=0.83mo  [11]  participate  in  the  transport. 

The  conductivity  a  and  Hall  coefficient  Rn  are  given  by 


o  =  e  (/ir«r  +  +  Hhh) 


(10) 


-rrH rnr  ~  rgijni  + 
e  (/ir^r  +  +  ^hh)2 


(11) 


where  rp,  r/y,  r^,  /ip,  fii,  /i/,,  «ri  and  h  are  T  and  L  electrons  and  hole  Hall  factors, 
mobilities  and  concentrations,  respectively. 

The  longitudinal  polar  optical  and  acoustical  modes  are  expected  as  a  dominant  electron 
and  hole  scattering  mechanisms.  The  Boltzmann  equation  is  solved  using  intraband  scat¬ 
tering  and  interband  scattering  between  T  and  four  L  minima  as  well.  Matrix  scattering 
elements  from  Nag  [14]  were  used  in  the  calculation. 

The  following  values  of  material  parameters  were  used  in  transport  calculations:  dielectric 
constants  e0  =  10.3,  e0 0  =  7.1  and  energy  of  longitudinal  optical  phonons  huiO—0.2\3 
eV  [11].  The  least  amount  of  information  is  available  about  the  deformation  potential  of 
conduction  and  valence  bands  A Ee  and  A £/,  needed  to  calculate  scattering  on  acoustical 
phonons.  Nag  shows  for  A  Ee  =  4el/  [14]  and  A  Eh  can  be  only  roughly  estimated  by  a 
method  proposed  by  Kranzer  [32].  We  have  used  A E1/,  =  lOeK  in  the  calculations. 
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In  this  paper  the  final  Hall  mobility  fin  =  \Ru&\  and  Hall  concentration  tin  =  |l/e/?//| 
are  calculated  and  compared  with  the  experiment.  The  product  \Rna\  calculated  first  for  the 
intrinsic  semiconductor  together  with  experimental  points  is  shown  in  Fig.  3.  The  essential 
influence  of  the  L  valley  transport  to  /i#  is  evident.  A  critical  analysis  of  interpretation  of 
Smith  [6]  and  Chern  [7]  HT  concentration  data  was  given  by  Su  et  al.  [28],  but  mobility 
data  are  in  both  papers  correct  and  are  shown  in  Fig.  2.  More  detailed  information  about 
mobility  calculations  is  in  [57].  The  concentrations  n//  are  evident  from  Figs.  4,  8. 


A.  Defect  Modeling 

We  briefly  outline  the  statistics  of  the  defects  which  is  used  in  this  paper.  The  approach 
and  used  symbols  are  the  same  as  in  the  Berding’s  paper  [26]. 

The  first  task  to  be  solved  is  a  decision  about  dominant  native  point  defects  which  control 
the  defect  statistics.  There  was  a  wide  discussion  whether  Te  vacancies  (VTe)  or  Cd  intersti¬ 
tials  ( Cdi )  dominate  as  donor  defects.  On  the  contrary  measurements  of  lattice  parameters 
exist  [22]  indicating  that  Vre  is  the  dominating  donor.  We  believe,  that  this  deduction  is 
not  fully  conclusive,  because  the  samples  on  which  the  measurement  was  performed  were 
prepared  by  quenching  from  the  three  phase  line  at  the  Cd  site  of  the  P-T  diagram.  It 
is  known,  that  such  samples  are  not  well  defined  (non  ideal  casting)  and  must  contain  a 
relatively  large  amount  of  Cd  precipitates,  which  can  influence  the  lattice  parameters  more 
strongly  than  Vre ■  The  ab  initio  calculations  [1]  prefer  Cdj  over  Vre.  The  acceptor  defects. 
Cd  vacancies  ( VCd )  are  usually  reported  and  the  calculation  supports  such  a  result,  too.  The 
donor- type  Te  antisite  ( Teed )  was  proposed  in  [1].  Its  influence  is  expected  to  be  strong 
at  Te-saturated  conditions.  Existence  of  two  P-N  lines  on  CT-C1  doped  samples  [31]  was 
interpreted  by  this  defect. 

In  the  first  approach  the  densities  of  neutral  native  defects  mentioned  above  are  given 
by  reaction  constants  [1] 


[Cdf]  = 


n, 


o'JfPc 


Cd 


I< 


exp) 


(E(cd,) + u(cd,rb)  +  s(cd,rb 


kbT 


(12) 
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TPh 

n,'/^P-  exp  (E(V4v)  +  U(VT')m)  +  S(VT,) 


vib 
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(15) 


E  is  relevant  defect  formation  energy,  energy  Umb  and  entropy  Svtb  determine  contributions 
to  the  vibrational  free  energy. 


K  =  kbT 3 


rncdhY' 2 

2'kTi2  ) 


n0  =  1.48  x  1022  cnfJ  is  the  number  of  unit  cells  per  volume,  kB  is  the  Boltzmann  constant 
and  m cd  is  the  mass  of  Cd  atom.  The  label  x  is  used  for  neutral  defect  overall. 

The  densities  of  multiply  ionized  defects  are  calculated  for  acceptors  and  donors 


[xz[ 

— 

Xx]pxi-exp 

[xz‘' 

— 

xx]  gX‘-  exP 

j-(znF-Ela-...-Exa) 
JT. j<(Ed  +  •••  +  Ea  -  zHf) 


(16) 

(17) 


where  Ea  and  Ed  are  the  acceptor  and  donor  one-electron  ionization  energies  and  jip  is  the 
Fermi  energy.  gXz'  {9x*‘)  is  the  degeneracy  of  the  defect.  The  index  z  means  the  ionization 
degree,  the  prime  and  bullet  correspond  to  negative  and  positive  species.  Experimentally 
determined  ionization  energies  of  native  defects  are  shown  in  Tab.l.  Important  are  especially 
deep  donor  levels  for  \rcd  anfl  acceptor  levels  for  Vre-  These  results  were  obtained  from 
EPR  measurements,  one  of  the  most  precise  and  reliable  methods  for  a  determination  of 
activation  energies  of  defects.  The  values  of  energies  and  entropies  established  by  theoretical 
calculations  for  defects  are  given  in  the  Tab. 2.  Based  on  the  given  forms,  the  ab  initio 
calculations  describe  the  properties  of  native  defects  in  CT  qualitatively  [1],  The  relations 
(16), (17)  present  a  simplified  statistical  model.  The  exact  relations  (see  e.g.  [33])  yield 
only  slightly  different  results.  The  deviations  from  the  experimental  data  are  however  too 
extensive  to  describe  the  relevant  quantities  in  a  practically  applicable  form.  In  this  respect 
we  suggest  a  variant  (semiempirical)  approach  in  this  paper  to  solve  the  task  of  native  defects 
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in  CT.  We  fit  the  experimental  high-temperature  a  and  Rh  optimizing  key  parameters  of  the 
ah  initio  output  (formation  energy,  vibration  free  energy)  by  parameters  to  hit  the  native 
defects  properties  in  CT  more  satisfactorily.  Henceforward  we  use  donor-type  (Cd,-  and 
T ecd)  and  one  acceptor-type  native  defect  (Vcd).  The  donors  are  described  by  (12)  and 
(14)  and  the  acceptor  is  described  by  (13).  The  set  of  native  defect  ionization  energies  and 
degeneracy  factors  which  are  used  in  the  calculations  is  collected  in  Table  1.  The  Fermi 
energy  is  obtained  solving  electric  neutrality  condition  in  which  T  and  L  electrons,  holes 
and  singly  and  doubly  ionized  native  defects  are  included. 

nL  +  nr  +  [V^d]  +  2  [l&]  =  p  +  [Cdj]  +  2  [Cd;']  +  [T eCd]  +  2  [Te£d]  (18) 

Comparison  of  concentrations  of  defects  calculated  with  the  help  of  eqs.  (12,13,14,15) 
shows,  that  Cd  vacancies  Vcd  play  the  dominant  role  as  acceptors  and  Cd  interstitials  as 
donors.  The  tellurium  antisites  Tec<*  can  have  a  smaller  influence  as  donors.  The  impact  of 
Te  vacancies  is  negligible. 


B.  Results  and  discussion 

The  same  band  parameters  from  III. A  which  fit  the  Hall  mobility  pn  =  \Rno\  were 
consecutively  used  for  the  modeling  of  native  defects.  The  electric  neutrality  condition  and 
the  Boltzmann  equation  were  solved.  The  theoretical  and  experimental  value  of  n«  =11/67?// 
are  plotted  in  Fig.  4.  It  is  necessary  to  bear  in  mind  that  nfI  in  the  case  of  occupied  L-band 
does  not  express  the  electron  concentration.  Measurements  performed  under  Cd  pressure 
(full  points)  exhibit  a  good  agreement  with  the  experiment.  On  the  contrary  the  discrepancy 
between  theoretical  values  and  experimental  data  of  measurements  done  under  Te  pressure 
(open  points)  is  slightly  larger  probably  due  to  a  higher  influence  of  Te  on  a  quality  of 
Mo  welded  contacts.  For  all  experimental  points  R//  <  0.  Corresponding  experimental 
and  theoretical  values  of  conductivity  a  in  dependence  on  Cd  or  Te  pressures  at  different 
temperatures  are  shown  in  Fig.  5. 
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Fig. 4.  The  experimental  (upper  triangles  -  GOO°C,  squares  -  700°C,  circles  800°C, 
bottom  triangles  -  900°C  and  squares  -  1000°C)  and  theoretical  Hall  concentrations.  The 
full  lines  show  the  fit  within  presented  model.  The  dotted  lines  are  results  when  the  Teed 
was  included  in  the  model.  The  dash  line  outlines  the  position  of  the  stoichiometric  (P  N) 
line  which  was  obtained  from  the  theoretical  calculations  based  on  the  presented  model. 
The  dash-dot  line  shows  the  position  of  the  line  of  congruent  sublimation.  The  arrows  show 
the  Cd  saturation  limit. 
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Fig. 5.  The  experimental  (upper  triangles  -  600°C,  squares  -  700°C,  circles  -  800°C, 
bottom  triangles  900°C  and  squares  -  1000°C)  and  theoretical  conductivity.  The  full  lines 
show  the  fit  within  presented  model.  The  dotted  lines  are  results  when  the  Teed  was  included 
in  the  model.  The  dash  line  outlines  the  position  of  the  stoichiometric  (P-N)  line  which  was 
obtained  from  the  theoretical  calculations  based  on  the  presented  model.  The  dash-dot  line 
shows  the  position  of  the  line  of  congruent  sublimation.  The  arrows  show  the  Cd  saturation 
limit. 
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The  fitted  donor  and  acceptor  energies  and  entropies  together  with  the  ab  initio  results 
[1,2G]  are  shown  in  Table  2.  We  present  total  formation+vibration  energies  only  because  the 
E  and  Umb  cannot  be  distinguished  in  the  fit.  The  agreement  of  theoretical  and  experimental 
values  is  very  good. 

The  effect  of  Te  antisite  Teed  cannot  be  estimated  within  our  model  definitely.  From 
the  point  of  view  of  the  HT  properties  at  the  Te  saturate  conditions,  where  Teed  should 
play  the  most  significant  role,  we  observe  a  strong  influence  of  the  hole  transport,  which 
results  in  the  decrease  of  the  Hall  mobility  in  this  region.  Due  to  the  low  mobility,  however, 
the  holes  do  not  revert  the  sign  of  Rh  which  stays  negative.  The  donor  Tea  compensates 
acceptor  Va  keeping  the  electron  high-mobility  transport  dominating  at  the  Te  saturation 
as  well.  Such  a  behaviour  is  evidently  not  observed  within  all  temperature  interval  where  the 
measurements  at  the  Te  saturated  side  were  performed.  We  can  thus  deduce  that  Teed  is  a 
minor  defect  above  600°C  and  does  not  manifest  in  the  electrical  measurements.  The  real 
concentration  of  charged  Tea  is  less  than  some  limit  value  which  prevents  Tea  to  prevail. 
We  have  estimated  such  a  value  within  the  parameters  in  Tables  1.  a  2.  by  a  value  E  >6 
eV  which  is  to  be  substituted  in  the  Eq.  (14).  The  open  question  remains  the  concentration 
of  the  neutral  Tea ■  The  ionization  energy  used  in  our  calculation  is  not  set  safely  and  the 
model  of  Teed  as  a  deep  donor  can  be  valid  as  well.  In  such  a  case  the  concentration  of 
neutral  Tea  could  increase  significantly  without  the  effect  to  the  electrical  measurements 
but  having  an  impression  to  the  solidus  curve  in  the  T-x  diagram. 

The  essential  effect  of  the  four  L  minima  in  the  conduction  band  to  the  native  defect 
statistics  is  outlined  in  Figs.  6,  7,  8,  9,  10  at  800°C.  Due  to  the  occupation  of  the  L  bands 
the  relevant  electron  concentration  is  not  expressed  satisfactorily  by  \\/eRn\.  If  higher 
minima  in  the  band  structure  are  occupied,  the  simultaneous  analysis  of  both  Ru  and  a 
describes  the  carrier  concentrations  and  mobilities  correctly.  Consequently  the  concentration 
of  native  defects  can  be  twice  or  more  higher  than  electrons  being  assumed  only  in  the  T- 
point.  The  dash  line  shows  the  concentration  of  electrons  and  holes  at  room  temperature. 
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Fig-6.  The  neutral  (x)  and  ionized  (’  ;  *)  native  defect  concentrations  together  with 
total  electron  (e’)  and  hole  (h*)  concentration,  all  at  600°C,  are  plotted  by  full  lines.  The 
experimental  (squares)  and  theoretical  (dotted)  \\feRH\  identical  with  that  in  Fig.  4  are 
involved.  The  dash  lines  show  room  temperature  e’  and  h*. 
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Fig. 7.  The  neutral  (x)  and  ionized  (’  ;  *)  native  defect  concentrations  together  with 
total  electron  (e’)  and  hole  (h*)  concentration,  all  at  700°C,  are  plotted  by  full  lines.  The 
experimental  (squares)  and  theoretical  (dotted)  \1/cRh\  identical  with  that  in  Fig.  4  are 
involved.  The  dash  lines  show  room  temperature  e’  and  h*. 
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Fig. 8.  The  neutral  (x)  and  ionized  (’  ;  *)  native  defect  concentrations  together  with 
total  electron  (e’)  and  hole  (h*)  concentration,  all  at  800°C,  are  plotted  by  full  lines.  The 
experimental  (squares)  and  theoretical  (dotted)  \\/cRH\  identical  with  that  in  Fig.  4  are 
involved.  The  dash  lines  show  room  temperature  e'  and  h*. 
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Fig.9.  The  neutral  (x)  and  ionized  (’  ;  *)  native  defect  concentrations  together  with 
total  electron  (e’)  and  hole  (h*)  concentration,  all  at  900°C,  are  plotted  by  full  lines.  The 
experimental  (squares)  and  theoretical  (dotted)  \l/eRH\  identical  with  that  in  Fig.  4  are 
involved.  The  dash  lines  show  room  temperature  e’  and  h*. 
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Fig. 10.  The  neutral  (x)  and  ionized  (’  ;  *)  native  defect  concentrations  together  with 
total  electron  (e’)  and  hole  (h*)  concentration,  all  at  1000°C,  are  plotted  by  full  lines.  The 
experimental  (squares)  and  theoretical  (dotted)  \l/eRu\  identical  with  that  in  Fig.  4  are 
involved.  The  dash  lines  show  room  temperature  e’  and  h*. 
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The  knowledge  of  equilibrium  concentrations  of  defects  at  HT  allows  to  establish  the 
crystal  properties  at  300K,  which  are  essential  for  technological  applications,  when  usually 
crystals  with  a  minimum  deviation  from  stoichiometry  are  demanded.  If  the  concentrations 
of  defects  at  HT  are  known  and  we  suppose  an  ideal  quenching,  we  can  determine  by  a 
solution  of  electric  neutrality  condition  at  300K  isoconcentration  lines  for  electrons  ni  and 
for  holes  h  (n*,  is  at  300K  negligible)  in  dependence  on  annealing  conditions  (T,  pCd ),  from 
which  the  crystal  was  quenched.  Isoconcentration  lines  are  shown  in  the  P-T  diagram  in 
Fig.  11.  The  region  of  pressures  p cd,  at  which  the  Fermi  level  of  the  sample  at  the  room 
temperature  is  at  the  same  position  as  it  is  in  the  intrinsic  semiconductor  can  be  seen  in 
this  plot.  Corresponding  line  (the  so  called  P-N  line),  for  which  the  following  expression 
was  found 

Pcd{P  ~N)=  4.5  x  105 exp(— 1.73  x  10 4/T){atm)  (19) 

is  given  in  Fig.  1.  Here  also  the  lines  corresponding  to  stoichiometric  composition  of  gas 
pc  and  liquid  pL  are  plotted.  This  P-N  line  roughly  gives  conditions  for  preparation  of 
samples  with  a  minimum  deviation  from  stoichiometry.  We  verified  its  position  by  aii- 
nealing/quenching  experiments  at  600°C,  700°C',  800°C,  900°C  under  different  Cd  or  Te 
pressures,  when  annealing  time  was  chosen  high  enough  to  establish  equilibrium  concentra¬ 
tions  of  defects  by  the  diffusion  process.  The  obtained  results  of  p cd  are  slightly  higher  than 
those  obtained  from  evaluation  of  HT  conductivity  and  Hall  effect  measurements  and  are 
lying  within  the  gray  area  marked  in  Fig.  12. 

From  the  point  of  view  of  preparation  of  single  crystals  with  a  minimum  deviation  from 
stoichiometry  it  is  necessary  to  know  the  curve  of  solidus  in  the  T-x  diagram.  This  can 
be  evaluated  from  equilibrium  concentrations  of  electrically  active  defects  established  from 
Hall  measurements  on  the  three  phase  line  [0,21,22]  or  directly  from  thermodynamic  mea¬ 
surements  of  Cd  or  Te  pressure  on  the  three  phase  line  [37,38,43].  Thermodynamic  and 
galvanomagnetic  measurement  of  Cd-rich  solidus  are  not  contradictory.  A  discrepancy  ex¬ 
isted  for  the  Te-rich  solidus, 
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Fig.ll.  F-T  diagram  with  isoconcentration  lines  for  electrons  np  and  holes  h 
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Fig.  12.  Partial  pressure  of  Cd  along  the  three-phase  loop  of  CdTe(S).  The  full  lines 
ps  PL ,  and  PG  represent  the  partial  pressure  over  stoichiometric  CdTe(S),  stoichiometric 
CdTe(L).  and  congruently  subliming  CdTe(S)  or  evaporating  CdTe(L),  respectively.  The 
Ps  was  found  based  on  results  of  the  presented  model.  The  grey  area  gives  the  region  of 
approximate  position  of  this  line  as  observed  from  our  annealing  and  quenching  experiments. 
P°  is  Cd  partial  pressure  over  pure  Cd. 
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especially  between  measurements  of  R//  [22]  and  thermodynamic  measurements  [37,45]. 
Based  on  results  of  our  experiments  it  is  possible  to  establish  the  solidus  curve  and  compare 
the  obtained  results  with  these  experiments.  The  arrows  in  Fig.  4  indicate  pressures  corre¬ 
sponding  to  the  three  phase  boundary.  The  solidus  curve  constructed  for  these  points  and 
for  the  whole  three  phase  line  is  shown  in  Fig.  13.  The  agreement  of  our  measurements  and 
those  of  Refs.  [37,43,45]  both  for  Cd  and  Te-rich  sides  is  good.  Comparison  our  results  with 
the  results  of  Wienecke  [22],  who  measured  R//  and  a  in-situ  in  dependence  of  temperature, 
Cd  and  Te  pressure  corresponding  to  the  three  phase  line  yields  a  good  agreement  on  the 
Cd  rich  side.  On  the  contrary  a  significant  discrepancy  can  be  seen  on  Te-rich  samples.  The 
interpretation  of  measurements  of  R#  and  o  on  the  Te  side  of  the  P-T  diagram  is  incorrect 
from  the  following  reasons:  The  assumption,  that  R u  is  not  influenced  by  electrons,  if  their 
concentration  is  roughly  one  order  of  magnitude  lower  than  intrinsic  concentration,  is  false. 
Even  due  to  the  fact,  that  the  value  of  mobility  ratio  b —(fie/nh)  ~  10  in  the  discussed  tem¬ 
perature  interval,  inversion  of  sign  of  R//  should  occur  at  580°C.  The  quantity  (p-n)  cannot 
be  found  from  measurements  in  weak  magnetic  fields  without  a  knowledge  of  mobilities  of 
electrons  /ie  and  holes  /!/,,  which  are  not  given  in  [22].  The  mobility  values,  which  strongly 
influence  the  obtained  concentrations  of  defects  are  not  presented  in  ref.  [23]  as  well. 

Further  contradictions  are  apparent.  The  authors  supposed,  that  the  electric  neutrality 
condition  p-n=N^  is  valid  in  the  whole  interval  of  temperatures.  This  means,  that  in  this 
temperature  region  acceptors  Na  are  fully  ionized.  This  can  be  reached  only  for  shallow 
acceptors.  In  such  a  case,  however,  the  slope  of  temperature  dependence  of  concentration 
(p-n)  determines  the  energy  of  formation  of  Vcd-  Using  the  data  of  Ref.  [22]  yields  the  value 
740meV,  which  is  substantially  lower  than  the  result  of  ab  initio  calculations. 

High-temperature  in  situ  Hall  effect  and  conductivity  measurements  are  presented  in 
this  paper.  The  results  show,  that  at  high  temperatures  two  types  of  electrons  participate 
in  the  transport:  Electrons  in  the  T  point  with  mr=  0.096m0  and  electrons  in  the  L  point 
with  mi  =  0.35m0.  New  experimental  values  of  energy  and  entropy  terms  of  native  point 
defects  in  CdTe  are  presented  and  comparison  with  the  ab  initio  calculations  show  a  good 
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agreement.  The  position  of  isoconcentration  lines  for  electrons  and  holes,  P-N  line  in  the 
P-T  diagram  and  the  solidus  line  in  the  T-x  diagram,  which  describe  well  the  experimental 
data  were  established.  The  dominant  preference  of  the  reported  research  is  comprised  in  the 
integral  description  of  various  features  of  CdTe  within  the  only  one  model.  The  most  of  the 
parameters  which  are  used  were  obtained  from  the  independent  experiments  and  the  fit  of 
several  of  them  do  not  deviates  significantly  from  their  theoretical  predictions. 


Fig.  13.  Tellurium  atom  concentration  in  CdTe  for  Cd  and  Te  -  saturated  conditions. 
Circles  represent  experimental  data  from  Ref.  [38]  and  squares  data  from  ref.  [37].  The  dash 
line  plots  the  gasus  stoichiometry  conditions. 
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V.  GALVANOMAGNETIC  PROPERTIES  BELOW  AND  ABOVE  THE  MELTING 

POINT 

To  prepare  high  quality  single  crystalline  CdTe  from  the  melt,  it  is  necessary  to  have 
sufficient  information  on  the  S-L-G  phase  equilibrium  and  on  the  L-G  and  S-G  equilibria. 
The  phase  equilibrium  S-L-G  was  studied  experimentally  both  on  the  Cd  side  [15,42  44] 
and  the  Te  side  [44,45]  and  theoretically  [1].  The  phase  equilibrium  S-G  is  relatively  well 
known  from  conductivity  and  Hall  effect  measurements  in  dependence  on  temperature  and 
pressure  of  one  of  the  components,  which  were  mainly  used  to  study  the  S-G  equilibrium 
experimentally  [15,0,7,22].  Due  to  the  fact,  that  multivalley  transport  of  electrons  has  to  be 
taken  into  account  (see  section  III),  a  new  interpretation  of  Hall  effect  and  conductivity  ex¬ 
perimental  data  received  so  far  was  performed(Section  IV).  The  least  amount  of  information 
is  available  about  the  L-G  equilibrium  above  the  melting  point.  The  equilibrium  pressure 
of  Cd  gas  over  the  liquid  characterizing  the  stoichiometric  composition  of  the  melt  is  shown 
in  the  p-T  diagram  in  Ref.  [46].  Usually  three  species  are  assumed  in  liquid  phase:  Cd 
and  Te  atoms  and  CdTe  molecules.  This  assumption  was  not  proved  experimentally.  Our 
growth  experiments  performed  by  gradual  cooling  of  the  melt  in  the  temperature  gradient 
show,  that  the  quality  of  the  crystal  can  be  closely  connected  to  the  history  of  the  CdTe 
melt.  Therefore  the  fabrication  of  high  quality  CdTe  crystals  represents  a  very  complex 
technological  problem  in  comparison  to  crystals  of  elementary  or  III-V  semiconductors. 

The  study  of  the  CdTe  liquid  is  interesting  also  because  it  embodies  an  extremely  dif¬ 
ferent  behaviour  which  differs  from  the  liquid  of  group  IV  and  III-V  semiconductors.  The 
basic  difference  reposes  on  the  fact,  that  strong  cohesive  forces  remain  after  the  destruction 
of  long-range  order  and  that  a  liquid  model  with  weak  interactions  is  inappropriate  in  CdTe 
in  comparison  to  e.g.  GaAs  [47].  Measurements  of  neutron  diffraction  [48]  confirm  the  exis¬ 
tence  of  the  short-range  ordered  network.  This  fact  is  confirmed  also  by  other  experiments 
performed  in  liquid  CdTe.  It  was  observed  [49,54],  that  the  dependence  of  the  electrical 
conductivity  on  temperature  in  the  liquid  CdTe  is  the  same  like  for  semiconductors  in  the 
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condensed  phase,  when  the  electrical  conductivity  rises  with  the  increasing  temperature. 
Such  a  behaviour  was  not  observed  experimentally  in  any  semiconductor  from  group  IV  or 
III-V,  where  the  electrical  conductivity  has  a  metallic  character  and  the  change  between 
the  semiconducting  and  metallic  behaviour  has  a  step-like  form.  Numerical  calculations  for 
liquid  CdTe  [2]  show,  that  metallic  behaviour  is  not  typical  for  CdTe,  probably  due  to  the 
iouicity.  A  higher  ionicity  of  Cd-Te  coupling  limits  the  delocalization  of  electrons,  which  can 
contribute  to  the  semiconducting  behaviour  of  the  liquid  CdTe. 

The  measurements  of  GM  properties  below  and  above  the  melting  point  were  performed 
in  two  stages.  In  the  first  stage  we  concentrated  on  conductivity  measurements.  Then  the 
measuring  cell  was  modified  and  improved  in  order  to  enable  Hall  effect  measurements. 

A.  Conductivity 

Our  DC  measurement  of  electrical  conductivity  was  performed  in  a  quartz  cell  in  the 
van  der  Pauw  configuration  and  spectrally  pure  graphite  was  used  to  prepare  contacts.  The 
temperature  was  measured  by  a  Pt-PtRh  thermocouple.  The  cell  was  connected  with  a 
CdTe  stack  and  both  were  filled  by  CdTe  powder.  After  melting  a  free  volume  in  the  stack 
was  approx,  the  same  as  a  volume  of  the  cell.  Due  to  the  fact,  that  it  was  not  possible 
to  control  the  Pcd  overpressure  above  the  sample,  both  the  solid  and  liquid  were  weakly 
depleted  of  Cd.  The  depletion  of  Cd  is  connected  with  a  change  of  the  concentration  of 
point  defects  (probably  Vcd)-  The  temperature  dependence  of  the  conductivity  was  measured 
after  melting  of  CdTe  up  to  1200°C.  Subsequently  the  cooling  was  performed  with  the  rate 
15°C/hour  at  temperatures  1200-1000°C  and  with  rate  30°C/hour  at  temperatures  1000- 
400°C.  This  process  was  followed  by  another  heating/cooling  cycle  between  400  and  1200°C. 
No  significant  difference  of  conductivity  between  both  heating/cooling  dependencies  was 
observed.  The  results  of  our  measurement  together  with  those  from  Ref.  [49,54]  are  shown 
in  Fig.  14.  Two  clear  slopes  can  be  seen  on  our  line  -  0.7  eV  in  the  solid  and  4.6  eV  in  the 
liquid. 
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The  equilibrium  concentrations  of  defects  in  the  whole  region  of  stability  of  solid  phase 
depend  on  the  temperature  and  gas  pressure  of  one  of  the  components.  It  would  be  optimum 
to  study  the  electrical  conductivity  on  a  sample,  at  which  the  position  of  the  Feriui  level  is 
close  to  the  intrinsic  Fermi  level,  it  means  close  to  the  solidus-stoichiometry-line  marked  in 
the  p-T  diagram  (Fig.  15)  as  Ps.  This  line  is  sometimes  called  the  P-N  line.  Due  to  the 
fact,  that  our  measurements  were  performed  in  a  wide  range  of  temperatures  it  was  possible 
to  compare  data  obtained  at  lower  temperatures  on  solid  CdTe  with  our  independent  in-situ 
galvanomagnetic  measurements  which  were  performed  at  defined  T  and  Pa  [57].  Based  on 
this  comparison  we  could  also  complete  the  Pa  into  presented  data  shown  as  the  points  in 
Fig.  15.  The  dash  line  plots  a  seeming  path  of  the  sample  through  the  P-T  diagram,  along 
which  the  measurements  were  performed.  We  notice  that  the  limit  Pa  — >  0.0015  atm  at  low 
T  is  not  a  real  effect  but  it  results  from  the  foreign  impurity  conduction  with  a  concentration 
5  x  1015cm-3  in  this  sample.  The  real  Pcd  ought  to  follow  a  line  nearly  parallel  with  the 
gas-stoichiometry-line  PG  as  it  is  well  expressed  at  1000/T  <  0.9A'-1  showing  a  depletion 
of  Cd  in  the  CdTe  solid  sample.  The  moderate  shift  of  the  points  away  the  P°  at  high 
temperature  is  probably  due  to  extrapolation  error  near  P° .  Because  of  a  fast  sublimation 
of  the  sample  near  P°  no  in-situ  galvanomagnetic  data  at  defined  Pcd  are  available  there  and 
the  comparison  with  the  points  near  Ps  was  done.  The  continuation  of  the  experimental  Pa 
into  the  liquid  was  calculated  via  solution  of  balance-mass-equations  where  stoichiometry 
liquidus  pressure  line  PL  was  applied. 
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Fig.  14.  The  conductivity  of  CdTe  vs  1000/T.  Our  results  are  shown  for  cooling  (full 
triangles)  and  heating  (open  circles),  respectively.  The  measurements  of  [49]  and  [54]  are 
included  as  (stars)  and  (open  squares),  respectively.  The  full  line  shows  theoretical  intrinsic 
electrical  conductivity  of  solid  Oi  according  to  [57].  The  labels  in  the  inset  point  to  the 
solid-gas  (S+G)  and  the  liquid-gas  (L+G)  phase  equilibrium,  the  medium  area  (S+L+G) 
corresponds  to  the  three-phase  S-L-G  equilibrium. 
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Fig. 15.  Partial  pressure  of  Cd  along  the  three-phase  loop  of  CdTe(S).  The  full  lines 
Ps,  PL ,  and  PG  represent  the  partial  pressure  over  stoichiometric  CdTe(S),  stoichiometric 
CdTe(L),  and  congruently  subliming  CdTe(S)  or  evaporating  CdTe(L),  respectively.  P°  is 
Cd  partial  pressure  over  pure  Cd.  The  full  diamonds  and  dash  line  plots  the  experimental 
path  through  the  region,  see  details  in  text.  The  labels  correspond  to  the  previous  figure. 
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The  sample  starts  to  melt  at  T  =  10G9°C  and  Pa  =  0.23  atm  and  is  completely  molten 
at  T  =  1090°C  and  Pa  =  1-04  atm.  As  it  can  be  seen  from  Fig.  15,  the  phase  transition 
solid-liquid  starts  in  our  case  close  to  the  intersection  of  the  PG  line  with  the  three  phase 
boundary.  Our  experimental  value  of  the  electrical  conductivity  at  the  melting  point  is 
«  10  Q-1cm-1,  the  measurements  from  Ref.  [49,54]  give  60  Q-1cm-1  and  90 
respectively.  These  differences  can  be  explained  by  the  fact,  that  the  measurements  [49,54] 
were  probably  performed  on  sample  with  a  larger  deviation  from  stoichiometry,  what  can 
be  deduced  from  the  lower  temperature  of  the  melting  point  in  their  measurements.  The 
conductivity  below  and  above  the  melting  point  of  CdTe  were  also  performed  by  non-contact 
EDDY  current  measurement  Ref.  [51].  The  values  of  conductivity  in  the  liquid  are  in  good 
agreement  with  the  results  of  our  measurements.  The  abrupt  increase  of  conductivity  from 
solid  to  liquid  state  in  [51]  compared  to  our  results  can  be  caused  by  a  smaller  deviation  from 
stoichiometry  in  their  sample,  due  to  the  minimized  free  volume  in  thp  ampoule.  Transition 
from  liquid  to  solid  state  in  this  case  occurs  close  to  the  maximum  melting  point  at  1092°C. 

B.  Hall  effect 

The  construction  of  Hall  cell  for  measurement  of  Hall  effect  in  liquid  CdTe  in  the  Van 
der  Pauw  configuration  represents  a  difficult  technological  problem.  The  keystone  of  our 
construction  of  the  experimental  setup  is  a  quartz  cell  with  four  pipe  outlets  for  a  connection 
of  contacts.  It  came  out,  that  most  metals  react  with  liquid  CdTe  resulting  in  a  metal 
dissolution  and  subsequent  strong  doping  of  the  liquid.  It  appeared  that  only  graphite 
could  be  used  as  contacts.  We  used  spectrally  pure  graphite  with  purity  5N.  Both  DC 
and  AC  methods  were  used  to  measure  a  and  R//.  It  appeared,  that  the  latter  is  more 
suitable  especially  for  measurements  above  the  melting  point  due  to  reduced  noise.  These 
measurements  were  performed  in  a  DC  magnetic  field,  when  AC  current  with  frequency 
600Hz  was  passing  through  the  sample.  The  measuring  cell  connected  with  a  stack  was 
evacuated  and  completely  filled  with  liquid  CdTe  after  heating  above  the  melting  point 
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temperature.  Several  measurements  of  transitions  from  solid  to  liquid  and  back  with  both 
increasing  and  decreasing  temperature  were  done.  The  goal  was  to  measure  the  temperature 
dependence  not  only  around  the  melting  point,  but  also  in  a  wide  range  of  temperatures 
(600-1250°C)  in  order  to  compare  the  obtained  values  of  Hall  coefficient  R//  with  data 
from  measurements  on  crystalline  CT  at  temperatures  (600-1000°C).  The  measurement  was 
performed  in  a  one  zone  furnace,  when  the  temperature  was  measured  and  controlled  by  a 
Pt-PtRh(10%)  thermocouple.  The  cell  and  the  stack  had  almost  the  same  temperature.  Due 
to  the  complexity  of  construction  of  the  measuring  cell  it  was  not  possible  to  maintain  defined 
pressure  of  one  of  the  components  above  the  solid  and  liquid  during  the  measurements.  Some 
free  volume  existed  above  both  the  solid  and  liquid  sample,  in  which  CdTe  dissociated  to 
components  Cd  and  Te2  existed  in  gas  phase.  The  state  of  binary  CdTe  existing  in  two 
phases  (S-L  and  L-G)  is  according  to  the  Gibbs*  phase  rule  determined  by  temperature  and 
pressure  of  one  of  the  components  above  the  solid  or  liquid,  only  on  three  phase  line  (S-L- 
G)  by  temperature.  Concretely  R h  and  a  depend  on  temperature  and  p cd  pressure  in  the 
region  of  stability  of  the  solid  and  only  on  temperature  on  the  three  phase  line. 

In  our  case  the  more  volatile  component  -  Cd  starts  to  leave  the  material  until  the  Cd 
pressure  pcd  reaches  the  values  close  to  the  line  of  congruent  sublimation. 

The  temperature  dependence  of  l/(eR//)  is  shown  in  Fig.  16.  The  arrow  indicates  the 
value  corresponding  to  the  melting  point.  The  sign  of  R//  is  negative  both  below  and  above 
the  melting  point.  The  temperature  dependence  of  electrical  conductivity  and  Hall  mobility 
is  given  on  Figs.  17  and  18,  where  the  melting  point  is  again  indicated  by  an  arrow.  The 
measurement  in  the  solid  state  wras  analyzed  by  the  way  described  in  the  preceding  sections. 

The  supposed  conditions  of  experiment  under  which  the  measurement  from  low  temper¬ 
ature  (600°C)  to  the  melting  point  (1092°C)  or  back  were  performed  are  very  close 
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Fig.lG.  The  conductivity  below  and  above  the  melting  point  of  CdTe.  The  arrow  indi¬ 
cates  the  temperature  of  the  melting  point. 


Fig. 17.  The  1  f(eRn)  below  and  above  the  melting  point  of  CdTe.  The  arrow  indicates 
the  temperature  of  the  melting  point. 
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Fig.  18.  The  Hall  mobility  below  and  above  the  melting  point  of  CdTe.  The  arrow 
indicates  the  temperature  of  the  melting  point. 
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to  those  discussed  for  the  experiment,  when  only  conductivity  was  measured  (F'ig.  15). 


C.  Discussion 

To  elucidate  the  interesting  properties  of  the  liquid  CdTe,  we  suggest  a  model  explanation 
which  describes  the  observed  experimental  data.  It  needs,  however,  a  detailed  and  critical 
discussion  and  also  more  extended  experimental  basis  to  be  verified  satisfactorily. 

We  discern  three  temperature  intervals  of  the  existence  of  the  liquid  with  a  distinct 
thermodynamic  or  electric  behaviour.  The  first  interval  (1069°C  <  T  <  1090°C)  corresponds 
to  the  melting.  The  three  phases  solid-liquid-gas  exist  simultaneously.  The  PCd  strongly 
increases  and  the  sample  is  significantly  enriched  by  Te  from  ~  1018cm-3  at  the  lower  T 
to  ~  102°cm-3  at  the  upper  limit  T.  Simultaneously,  the  conductivity  is  not  substantially 
influenced  by  that.  This  means  that  Te  is  not  our  electrically  active  species  here. 

In  the  second  interval  (1090°C  <  T  <  1160°C)  the  CdTe  is  fully  molten.  The  escape 
of  Cd  from  the  sample  stops.  The  liquid  behaves  continuously  as  a  semiconductor  with  a 
highly  thermally  activated  free  electron  concentration.  We  interpret  the  strong  increase  of 
a  characterized  by  the  energy  4.6  eV  as  a  formation  of  electrically  active  Cd.  We  expect 
that  the  short-range  ordered  medium  coexists  with  randomly  distributed  Cd  and  Te  atoms. 
In  the  ordered  form,  the  covalent  binding  between  II- VI  compounds  survives  resulting  in 
the  semiconductor-like  behaviour.  The  disordered  atoms  can  be  regarded  as  defects  in  the 
lattice  being  both  donor-  and  acceptor-  like,  generally.  From  point  of  view  of  the  most 
frequently  used  electrically  active  native  defects  in  the  solid,  which  are  Cd  interstitial  Cd/ 
and  Cd  vacancy  Vcd,  we  recommend  to  describe  the  process  to  be  a  dissociation  of  the 
ordered  liquid  CdTe  according  to  the  reaction 

CdTe  — ^  TeVcd  T  Cd/,  (20) 

where  the  final  products  are  the  vacancy  in  the  correlated  ordered  lattice  Vcd  and  free 
unbound  Cd/.  The  observed  slope  in  a  and  Rh  is  in  a  good  agreement  with  first-principle 
calculation  of  this  energy  in  the  solid  which  comes  out  4.2  eV  [1]. 
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While  Cdj  is  the  electrically  active  divalent  donor,  the  Vcd.  in  the  ordered  soft  system 
relax  due  to  the  Jahn-Teller  distortion  of  the  surrounding  Te  atoms.  Such  a  relaxat  ion  shifts 
the  acceptor  level  above  Fermi  energy,  resulting  in  a  formation  of  neutral  Te-enriclied  region. 
Perspectively,  the  shift  up  could  convert  the  Te  level  to  become  donor. 

In  the  third  interval  (T  >  1160°C)  the  a  and  Rh  increase  moderately  in  comparison 
with  the  lower  T  resulting  to  the  slope  0.8  eV.  We  interpret  such  a  behaviour  as  a  hybridiza¬ 
tion  of  the  conduction  band  with  the  Cd  or  Te  impurity  band.  The  liquid  conies  to  be  a 
seinimetallic  degenerate  system.  The  defects  are  not  fully  ionized  and  the  system  cannot  be 
described  successfully  within  a  simple  semiconductor  model.  We  cannot  distinguish  from 
our  measurements  if  the  Cd  or  the  Te  impurity  band  is  responsible  for  such  effect.  In  view 
of  the  good  agreement  of  o  with  the  data  of  [49,54]  where  their  experimental  setup  resulted 
probably  into  increased  contents  of  Te  in  the  liquid,  we  prefer  hybridization  with  the  Te 
impurity  band. 

Evaluation  of  all  data  is  complicated  by  the  fact,  that  Cd  or  Te  pressure  in  the  measuring 
cell  are  not  well  defined  in  the  present  setup.  Improvement  in  precision  and  understanding 
of  results  of  GM  measurements  near  the  melting  point  of  CdTe  can  be  reached,  if  a  new 
setup  enabling  measurements  in  more  defined  conditions  is  constructed. 


VI.  CONCLUSION 

This  final  report  comprises  a  complex  view  of  defect  structure  of  CdTe  based  oil  high 
tenperature  GM  measurements  below  and  above  the  melting  point,  when  a  model  based  oil 
Cd  interstitials  and  Cd  vacancies  as  dominating  native  defects  is  applied.  The  possible  role 
of  Te  antisite  within  the  model  is  discussed. 
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TABLES 


TABLE  I.  The  ionization  energy  and  degeneracy  of  Cdj  and  Va ■  All  defects  are  divalent.  The 
two  energies  for  the  Cdi  correspond  to  the  two  tetrahedral  interstitial  sites  (ab  initio).  The  Vcd 
energies  are  referenced  in  the  Table,  other  data  are  acording  to  [1,39] 


Defect 

El  (meV) 

9i 

I52(meV) 

92 

Cdi 

0;  210 

2 

170;  360 

1 

Vcd 

50  [35] 

2 

470  [34] 

1/3 

Teed 

0 

2 

390 

1/3 

VTe 

1400  [40] 

2 

1800  [41] 

1 

TABLE  II.  The  ab  initio  calculated  energy  and  entropy  terms  and  the  respective  values  which 
were  used  in  the  fit  of  experimental  data.  The  underlined  parameters  were  optimized.  The  two 
energies  for  the  ab  initio  Cdj  correspond  to  the  two  tetrahedral  interstitial  sites. 


ab  initio 

fit 

Defect 

E  +  Umb(eV) 

S(h) 

E  +  Uvlb(cV) 

S(kb) 

Cd, 

0.96;  1.29 

8.8 

0.96;  1.29 

11.1 

VCd 

3.24 

-8.2 

3.35 

-5.6 

Teed 

4.72 

-17 

>6 

-17 
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